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Granular Segregation in Tapered Rotating Drums 
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We study the granular segregation in a tapered rotating drum by means of simulations. In this geometry, both 
radial and axial segregation appear, the strength of which depends on the filling fraction. We study the effect 
of the drum's tapering angle on the segregation speed and show that, coherently with several recent studies, the 
axial segregation is due to the combined effect of the radial segregation and the shape of the drum. We show that 
the axial segregation behaves in a way analogous to the one found in chute flows, and corresponds to previous 
theories for shallow gravity driven surface flows with diffusion. By means of this analogy, we show that tapered 
drums could be the simplest experimental set-up to obtain the free parameters of the theory. 

PACS numbers: 45.70.Mg 


It has been hinted that axial segregation in rotating drums 
is mainly the result of the shape of the free-surface, the thin 
layer of flowing particles in the drum ffl . The free-surface 
is not flat, but slightly curved along the axial direction, with 
peaks near the parallel walls of the drum. Due to the fric¬ 
tion with them, the particles on the sides are lifted more than 
the particles in the centre, thus creating a convective flow in 
the drum that induces the axial segregation. How to harvest 
this relation between flow and the segregation pattern remains 
an open question: a question whose answer could have far- 
reaching practical implications. 

Motivated by the results of - where axial segregation 
of binary mixtures in a tapered rotating drum is reported albeit 
not explained - and how they relate to our own results |j7[] 
- where the shape of the rotating drum is used to speed up 
and control the axial segregation - we study the segregation 
pattern in tapered rotating drums by means of simulations, and 
see to what extent it is possible to control it by modifying the 
drum’s shape. 

It has recently been shown that, contrary to the common be¬ 
lief that 3D rotating drums are the sum of independent quasi- 
2D slices interacting only via diffusion, the shape of the drum 
does influence the surface flow creating a slow axial drift 0]. 
This drift is caused by the axial slope of the drum’s walls, 
which determines the 3D structure of the flow. If a bi-disperse 
mixture is used, the particles will follow the flow structure and 
will segregate accordingly. 

In the particular case of a tapered drum, when mixed parti¬ 
cles enter the flowing layer in the upper part of the avalanche, 
gravity drives them towards the small side of the drum before 
returning them to the large one. Small particles can percolate 
trough the voids created by the big ones leaving the flowing 
layer before it changes its direction. This leads to a net trans¬ 
port of small particles to the small side of the drum and of 
large particles to the big side: axial segregation occurs. Since 
the shape of the flowing layer depends on the shape of the 
drum, one would also expect it to affect the segregation. In 
this paper we show that this is indeed the case and we study 
how it happens. 

The coupling of axial segregation with the velocity pro¬ 
file has been already reported for the classical rotating drum 
0 ], In the tapered rotating drum the same coupling is in play 


but with a fundamental difference: the shape of the drum in¬ 
creases the curvature of the flow lines on the flowing layer 
thus transporting more small particles than in a parallel drum, 
and in contrast to the latter, in a definite direction. 

The structure of the paper is as follow. First, the drum ge¬ 
ometry and and simulation details are explained. Second, we 
compare our results for half-filled drums with experiments to 
validate our model. Third, to minimize the effect of radial seg¬ 
regation, a lower filling fraction is used and total segregation 
is observed. The effect of the drum’s tapering angle on the 
speed of the segregation is also studied. Finally, comparisons 
will be drawn with simulations of mixtures in 2 and 3D chutes 
and with a theory for shallow gravity driven surface flows with 
diffusion. 
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FIG. 1: (Color online) Snapshots for half-filled system at different 
positions along the rotation axis (from left to right) and different 
times (from top to bottom, the number of revolutions indicated in 
the left). Small particles dark, large particles light. One can see that 
first the radial segregation appears, and this slowly and continuously 
transforms into axial segregation, that here is visible in the domi¬ 
nance of small particle in the small side of the drum. 
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SIMULATIONS 

Our simulations are in a slightly different geometry than 
that of [51 61 (shorter drum) to speed up the simulation time. 
We use an equal volume mixture of 5 and 10mm radii parti¬ 
cles in a tapered drum 150mm length, with caps of radii r = 50 
mm (at z = 0mm) and variable R (at z = 150mm), in the left 
and right side, respectively. Their relation define the angle of 
the drum, such that R = r(l +tan(a)). In the simulations, par¬ 
ticles’ contact forces are represented using a standard spring 
dash-pot model 0] as implemented in 0. The angular veloc¬ 
ity is set to 71 s , i.e. 30rpm, which implies an average Froude 
number of Fr = co 2 (r + R)/2g ^0.1. The system is either at 
50 or 25% filling fraction. As initial condition the particles 
are arranged in a square lattice, alternating one big particle 
and eight small to obtain an initial mixed state. Snapshots of 
the system for different times and axial positions are shown in 
Fig.Q](cf. Fig. 5 in Ref. 0 for an experimental comparison). 

The simulations were performed using Mercury-DPM [9[] 
with the same material properties used in our previous re¬ 
search 0, since they proved to reproduce the experimental 
results. 

Shape of the Free Surface 

In this section we show the three-dimensional structure of 
the flow. We focus on three different characteristics: the level 
of the surface along the axis of rotation (Fig. [3, the veloc¬ 
ity field on the free surface (Fig. [3} and the density profile, 
as seen by a vertical cut in the centre of the drum (Fig. ED. 
Together, they provide a rather complete way to visualise the 
structure. Data are averaged for 30 consecutive snapshots (one 
every 0.8s) and coarse-grained to obtain a continuous function 
following the procedure of Ref. 0. The weight function is 
a Lucy polynomial 0 with a width of 2d. Unless otherwise 
noted, every coarse-graining presented here is done with the 
same parameters. We measured at different times during the 
evolution of the system, and the results do not significantly 
change, so we conclude that what we observe depends only 
on the geometry and not in the degree of segregation. 

The free surface along the axial direction shows level dif¬ 
ferences. This phenomenon, to the best of our knowledge, 
was first described by Zik et al. for a drum modulated along 
the axial direction with a sinusoidal function 0. This can 
be seen in our simulations from Fig. |T] or looking at the iso¬ 
lines of the density at different axial positions, in Fig. [2] This 
level difference along the axis makes particles to move along 
the gradient of the surface, that is as particles go down the 
avalanche, they move forwards and backwards in the axis per¬ 
pendicular to the page. 

Particles enter the flowing layer with a small velocity and 
they are mainly directed towards the small side of the drum. 
This is due to the tapered geometry: when entering the flowing 
layer, the particles in the large side of the drum are higher 
than those in the small side of the drum, thus creating a flow 


towards the small side. This is inverted when particles reach 
the lower part of the avalanche: particles in the large side of 
the drum are lower than those in the small side. Increasing the 
angle of the drum, increases the curvature and the magnitude 
of this flow. This can be seen in Fig. [3 for three different 
angles. The surface is defined as the plane whose projection 
can be seen in Fig. [2] (solid line). 

When a mixture of particles enters the avalanche in the up¬ 
per part, they will segregate as they travel downslope. If the 
drum is large enough, some small particles will reach the core 
before the flow changes its direction, creating a net transport 
of small particles. (Large enough here means in relation to 
the segregation speed of the mixture. For size-driven segre¬ 
gation, the smaller the size ratio of the particles, the slower 
segregation is, and thus the drums needs to be larger.) The 
large particles, since they tend to stay on the surface, can flow 
back to the large side of the drum, following the flow lines. In 
this way, axial segregation occurs. 



x[m] 

FIG. 2: Contour lines for a density p* = p par ticles/ 10 at three dif¬ 
ferent axial positions for the same systems as in Fig. 1. The level 
difference can be seen mostly in the upper part of the avalanche. One 
should think of the curves as being displaced in the direction normal 
to the page. Compare the qualitative agreement with the model of 
Ref. [12]. Note that since the flow velocity increases along z, the 
blue curve shows a slightly more S-shaped profile. Solid line shows 
the plane where the velocity is measured for Fig.[3and dashed verti¬ 
cal line the plane where the density is measure for Fig. [4] 

Finally, the packing fraction of the system gives us an in¬ 
sight on the inner structure of the flow, and shows the depen¬ 
dence of the free surface on the drum’s angle. Fig. [4] Despite 
the presence of segregation, the density is uniform during the 
evolution of the system, with the densest region at the bottom 
of the drum. The density is mostly homogeneous along the ax¬ 
ial direction with increased spatial fluctuations as a increases. 

The density remains constants during the segregation pro¬ 
cess, and is maximum in the lower region of the drum for all 
the angles studied. We confirmed that the density does not 
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FIG. 3: Average velocity field of the particles’ free surface for the tapered drum (from left to right, tapering angle a = 10°, 15°, 20°), the colour 
representing the particles’ mean velocity in m/s. The length of the component of the vector field parallel to the axis of rotation is increased 
to stress the effect of the tapered geometry on the flow. As expected, the convective flows increases its amplitude and reach with the angle of 
tapering. 


vary by measuring it at different times during the segregation 
process. 


Radial and Axial Segregation 

To corroborate that the proposed mechanism is responsible 
of the axial segregation, we measure the radial and axial seg¬ 
regation in the drum as a function of time. In order to have 
order parameters for both, the radial segregation is defined as 
follows: 


^radial —^Ys<z* /Nts + Nl>z* / NtL 1 


( 1 ) 


where N s< ( L> y* is the number of small (large) particles whose 
radial coordinate is smaller (larger) than the half volume ra¬ 
dius r* (note that this radius depend on the axial coordinate), 
Nt s (l) being the total number of small (large) particles. That 
is, we count the fraction of particles that are in the expected 
side of the drum. When <yradial = 1 all the large particles are in 
the outside and the small in the core of the drum. It must be 
noted that, due to the flowing layer, the radial segregation can 
never reach the maximum value. To measure the axial segre¬ 
gation we use a quantity that also works with banding (which 
appeared in one of the systems studied) 111311 : 


W«) = W 

/ Jo Cs{z) 

where c s (z,t) is the number of small particles along the axis, 
and c s (z) the number of small beads that correspond to a slice 
of the drum at position z for the totally mixed state. It must 
be noted that since we are simulating a tapered rotating drum, 
the mean number of small particles changes with z since the 
cross section of the drum increases accordingly. 

Figure^shows the temporal evolution for each kind of seg¬ 
regation for three different angles. As expected, the radial 
segregation sets in quite fast, essentially after four turns, but it 


saturates soon to a roughly constant value. On the other hand, 
the axial segregation is slower but steadily increases beyond 
the point of saturation for the radial segregation. Eventually 
(around the 10th revolution) the axial segregation causes the 
value of q ra diai t° decrease slightly, since in the small side of 
the drum the number of large particles decreases greatly. 

If a is increased the segregation speed should increase 
since the level difference on the free surface is also increas¬ 
ing 0. This is indeed the case as |5l] shows experimentally 
(see in particular Fig. 8). The same principle can be seen in 
other configurations such as the spherical drum |14f|, or the 
double cone configuration 0. 

To better understand the role of geometry on the segregation 
process, several simulations were run for different a. It is 
expected that the radial segregation will not depend on the 
geometry of system. However, for the axial segregation, since 
it depends on the velocity profile and on the geometry of the 
flowing layer, one expects to find a strong dependence on the 
angle of the drum since this determines the level difference in 
the flow. Indeed this is the case as Fig.0shows. 


Less than Half Filled Drum 

Our previous result shed light on the relation between a 
and the axial segregation, but this is strong influenced by the 
radial segregation in half filled drum. So large in fact, that 
at this filling fraction some systems segregate the other way 
around, that is, large particles in the small side of the drum, 
as reported in 151] and confirmed by some of our simulations 
(data not shown). We ran simulations at 25% filling fraction, 
for different a, and we observe indeed that the axial segrega¬ 
tion is bolder. What is more, the radial segregation disappears 
almost completely since the mixing at this filling fraction is 
maximum 111611 , In this configuration, a pure phase of small 
particles appears in the small side of the drum and a steady 
state for the segregation is reached, while increasing the angle 
of the drum increases the segregation speed (see Fig. 
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FIG. 4: (Color online) Packing fraction on the middle plane of the drum for increasing drum angles a = 10°, 15°, 20° , from left to right. The 
increasing slope of the free surface can be seen from left to right, as a increases. 



FIG. 5: Axial (solid-lines) and radial (dashed-lines) segregation as 
a function of time for the tapered rotating drum. From light gray to 
black in increasing angle of tapering. As it can be seen, there is no 
correlation between the angle and the radial segregation. However, 
axial segregation shows a clear correlation with a. 

This behaviour is strikingly similar to that of bi-disperse 
chute flows, either in 2D [17] or 3D liTsil . Numerically, those 
are simulated with a box with periodic boundary conditions 
in the flowing plane, and a rough inclined bottom. Depend¬ 
ing on the inclination of the bottom with respect to gravity, 
a steady flow can be obtained. This is a straightforward way 
to study segregation for a fixed velocity profile, however their 
experimental realisation is rather complex. 

Figure[6]shows the displacement of the centre of mass of the 
large particles normalized by its asymptotic value, for three 
different drum angles 

\ Z C m{°°) — Z C m{t) ^ 

Az = -7—7- 7K\ ) I 3 ) 

where z cm (t) is the centre of mass position on the axial direc¬ 
tion for the large particles. As a increases, the system segre¬ 
gates faster since the shear rate increases with it. One can see 
that the drum’s angle influences the segregation in the same 


way as the size ratio between the particles lfl8[fl9l l or the vol¬ 
ume fraction between the two species 0. 



FIG. 6: Displacement of the centre of mass of large particles for 
systems with different a, normalized by the asymptotic value Az in 
a log-linear scale. Straight lines correspond to the exponential fit. 
Compare with Ref. d- If plotted in a linear scale, one obtains 
results as in Ref. I 1 18ll. 


This is valid for a wide range of drum angles. We define 
the segregation speed as the inverse of the exponent of the 
exponential fit Az = exp(— t/z). In general, the segregation 
speed is a monotonic function of the angle, as one can see in 
Fig. [7] the line being a linear eye-guide. A simple estimation 
tells us that the distance a particle travels on the free surface is 
L s oc tan(a). If the particles travel twice as much during one 
avalanche they will segregate twice as fast, and this would 
explain the linear behaviour of the segregation speed with the 
angle. 

Since we have this qualitative agreement between two com¬ 
pletely different systems, the obvious next step is to see how 
well the theoretical predictions developed for the chute flow 
work on the tapered drum. 
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FIG. 7: Normalised segregation speed measured as the inverse of the 
typical time, To/t obtained from the centre of mass displacement on 
Fig.[6]for different tapering angles (To is the typical time for a =2°). 
The solid line is a linear eye guide. 


Comparison with theory 


Let us look at the phenomenon from another point of view. 
Since the free surface is at an angle with respect to the hori¬ 
zontal (going up towards the large side of the drum, see Fig. 
2) the projection of gravity on the free surface is non-zero. 
This translates into an effective gravity towards the small side 
of the drum. From the theory for shallow granular avalanches 
|2(3, ;21] we know the time evolution for the concentration of 
each species, and we should be able to apply this result to the 
tapered drum. 


Whereas the chute flow develops vertical segregation as it 
moves downslope, the rotating drum develops axial segrega¬ 
tion as time passes. This is due to the structure of the flow 
as seen in Fig. Q] The tapered drum turns the kinetic siev¬ 
ing along the flow lines into axial segregation at the system’s 
level. According to this simple analogy, small particles should 
go to the small side of the drum, while large particle will go 
towards the large side. Furthermore, since the free-surface is 
a geometric characteristic of the flow and does not depend on 
the degree of segregation of the mixture, the tapered drum will 
segregate the particles into stable configurations in the same 
manner that a chute flow does. 


The temporal evolution of the concentration is something 
particularly difficult to measure experimentally and, to the 
best of our knowledge, there is only one realisation lf22tl . Our 
results show a strikingly good fit with the theory. Fig. [8] For 
the sake of space, we refer the reader to the original paper for 
the derivation of the solution 0 and to the appendix for the 
explicit formula we use. Note that for a tapered drum the half 
volume height is at 750 = L/y / (2), the spatial coordinate of 
the theoretical solution has been scaled to account for this. 


The results are remarkably similar concerning both the evo¬ 
lution of the concentration and the spatial distribution. How¬ 
ever, we cannot fit both the temporal evolution and the steady 
state solution. We chose to fit the Peclet number Pe (a dimen¬ 
sionless number that quantifies how large the segregation is 


compared to diffusive remixing) to reproduce the evolution of 
the system rather than the steady state. Does this mean that 
the Pe in the tapered drum changes with time? Answering this 
question goes beyond the scope of this paper. 

It must be noted that we are comparing the results of only 
one realisation, and our data points are not averaged on time, 
which would decrease the noise and give a better fit, but we 
prefer to show the raw results to highlight the quality of the 
data that is possible to obtain with this set-up. We confirmed 
that the same qualitative behaviour holds for different angles 
and drum lengths (data not shown). 

Finally, it must be noted that Gray & Chugonov’s theory 
can be used to explain the centre of mass movement reported 
in 1121 llSIl' together with our results. Consider the solution for 
the concentration of small particles in the chute, <p(x,t). To 
obtain the centre of mass of the small particles one integrates 
the concentration times the position, z C m(t) = Jq z 0(z,t) dz , 
which can be solved numerically in a very straightforward 
way. It is easy to verify that the predicted solution shows an 
exponential decay for small values of Pe ~ 4, but seems to 
show two different slopes for larger values of it. In this way, 
Pe can be measured not only for the steady state solution, as 
done in El but from the evolution of the system. The results 
presented here indicate that both may be different. 

A detailed study of Pe for different system parameters is be¬ 
yond the scope of this paper. Our objective being to show how 
well the theory describe the evolution of the system and how 
the tapered drum could become a trivial experimental set-up 
where to measure the Peclet number. A detailed study of this 
should be hopefully a subject to be addressed both experimen¬ 
tally and theoretically in the future. 


CONCLUSIONS 

We have studied the size-driven segregation of bi-disperse 
mixtures in tapered rotating drums. The relation between ve¬ 
locity field, free surface, and the segregation pattern was stud¬ 
ied. The results confirm that level differences on the free sur¬ 
face create axial flows, which in turn cause axial segregation. 
For the first time we report total axial segregation in this par¬ 
ticular geometry and we studied its dependence on different 
parameters; in particular, we showed that the axial segrega¬ 
tion speeds up with the angle of the drum. 

We have shown that tapered drums can be used to determine 
the Peclet number used in previous theories from a straight¬ 
forward experimental set-up. This has been done before but 
the experiments are very difficult to set up “with both well 
controlled initial conditions and a steady uniform bulk flow 
field” [220. And indeed, to the best of our knowledge, the 
results from Wiederseiner et al. are the only ones that give 
the Peclet number experimentally. On the contrary, to obtain 
well controlled results on a tapered drum is extremely easy, 
for example, using the same set-up that motivated this study 

j|. 

Besides these theoretical implications, the results presented 
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FIG. 8: Left panel, concentration of small particles along the normalized rotation axis as a function of time for our numerical experiments with 
a = 10° and L = 200mm. Right panel, theoretical solution using the results of Ref. (23tl with Pe = 14. 


here have direct applications in the separation and sorting of 
granulates in areas ranging from powder handling in the phar¬ 
maceutical industry to ore milling in the mining industry and 
shed light on the important role of geometry in segregation 
processes, not only for rotating drums but everywhere tridi¬ 
mensional flows appear. 
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Appendix I 

The general solution for the concentration of small particles 
can be written as 


0 


1 

2 


2 dco 


where co is decomposed onto a steady and a time-varying part 
co = (O s + CO T . Their analytical forms is given by 

_ v (P ,_sinh(^/2)_ _ sinh((|-4o)/2) 

' MSo) s . nh(( £ o/2) sinh ( 4 o/ 2 ) ’ 


rai(Pe 2 + 4 n 2 7 l 2 ) ^ ^ ^ 

In our solution, we take only the first hundred terms of the 
series. 
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